Comb-based spectroscopy with synchronous sampling for real-time averaging

ABSTRACT

A method of comb-based spectroscopy with synchronous sampling for real-time averaging includes measuring the full complex response of a sample in a configuration analogous to a dispersive Fourier transform spectrometer, infrared time domain spectrometer, or a multiheterodyne laser spectrometer. An alternate configuration of a comb-based spectrometer for rapid, high resolution, high accuracy measurements of an arbitrary cw waveform.

The present disclosure claims priority to U.S. Provisional Patent Application No. 61/243,554, filed Sep. 18, 2009.

BACKGROUND

The present disclosure relates to a comb-based spectroscopy.

Fourier transform spectroscopy (FTS or FTIR) has long been a workhorse system for research and industry. Recent advances in spectroscopy using dual frequency combs offer an interesting new approach to FTS. Frequency combs are stabilized pulsed lasers that produce a comb of well-defined frequency lines (hence the name frequency comb) in the frequency domain, equivalent to many cw lasers. They can also be considered in the time domain where they emit a train of optical pulses with a well-defined carrier frequency and repetition rate. Frequency comb sources are attractive for spectroscopy since they can be used, in principle, to combine the benefits of heterodyne laser spectroscopy with broad spectral coverage for high resolution, broadband measurements of samples. In a dual-comb spectroscopy approach, one comb probes a sample and a second comb acts as a local oscillator to read out the response. This dual-comb spectroscopy approach can be viewed as a form of infrared time-domain spectroscopy (TDS) analogous to THz TDS, or as a massively parallel multi-heterodyne laser spectrometer. It can be used to perform spectroscopy on a passive sample or on an active cw source.

SUMMARY

A method of comb-based spectroscopy with synchronous sampling for real-time averaging according to an exemplary aspect of the present disclosure includes generating a source comb that transmits a source pulse train through a sample to provide a signal pulse train; generating a local oscillator (LO) comb that transmits a LO pulse train at a different repetition rate than the source pulse train; establishing a base coherence between a first tooth of the source comb and the LO comb by phase locking each to a first CW reference laser and a second tooth of the source comb and the LO comb by phase locking each to a second CW reference laser with the number of intervening teeth for the source comb being N1 and the number of intervening teeth for the LO comb being N2 such that N2/(N1−N2) is an integer N, wherein each Nth pulse from the LO pulse train will exactly overlap each Nth+1 pulse from the source pulse train; detecting the signal pulse train through linear optical sampling against the LO pulse train in which consecutive samples of an overlap between a signal pulse and a LO pulse from the LO pulse train yields a measurement of the signal pulse; digitizing the overlap between the signal pulse train and the LO pulse train synchronously with the LO pulse train; and real time summing of every Nth digitized sample to generate an averaged signal pulse.

BRIEF DESCRIPTION OF THE DRAWINGS

Various features will become apparent to those skilled in the art from the following detailed description of the disclosed non-limiting embodiment. The drawings that accompany the detailed description can be briefly described as follows:

FIG. 1A is a simplified picture of a dual-comb spectrometer of spectroscopy of a sample. The pulse train from a source comb passes through the sample (and is possibly split as well to pass around the sample to yield a time-multiplexed normalization signal), then combined with a Local Oscillator (LO) pulse train. Depending on the configuration, a tunable bandpass filter is employed to satisfy Nyquist conditions. The heterodyne signal between the source and LO is detected and digitized to yield the complex gas sample response. In this configuration where the LO bypasses the sample, both the phase and magnitude of the sample response are measured;

FIG. 1B is a time-domain description of the coherent dual-comb spectrometer in which the source and LO emit pulse trains in time. Since they have different repetition rates, the relative source and LO pulse timing increments by ΔT=Δf_(r)/(f_(r,S)×f_(r,LO)) with each sequential pulse, where f_(r,S) is the source comb repetition rate, f_(r,LO) is the LO comb repetition rate, and Δf_(r)=f_(r,S)−f_(r,LO) is their difference.

FIG. 1C is a representation of actual digitized data of the photodetector voltage corresponding to the overlap of the source and LO pulses. Each data point corresponds to the product of the source and LO pulses, integrated over the pulse duration. The detected signal is plotted verses effective time by assigning the data point spacing to be exactly ΔT. Note that the data is equivalent to an interferogram, or cross-correlation, formed as pulses from the source and LO pass through each other. The large “centerburst” ranging from −3 ps to +3 ps corresponds to the overlap of the LO and source pulses. Also visible at 14 ps is a weak ringing containing absorption information from an HCN gas sample.

FIG. 1D is an equivalent frequency domain description of the coherent dual-comb spectrometer. In the frequency domain each comb creates an array of discrete teeth with slightly different spacings (repetition rates).

The right side of FIG. 1D represents the rf signature created by the two combs, each having slightly different repetition rates (tooth spacing), beating together to yield a third comb in the rf with spacing Δf_(r)=f_(r,S)−f_(r,LO). This rf comb is related to FIG. 1C simply by a Fourier transform a concatenated series of such interferograms and rescaling of the frequency axis. Note that for our experimental conditions the rf comb in would contain ˜4000 teeth;

FIG. 2A are example data traces illustrating the basic time domain signal and the ability of coherent signal averaging to improve signal-to-noise ratio (SNR). Here, the measured time-domain data (interferograms) are for a sample of 25 Torr of hydrogen cyanide (HCN) gas around 1539 nm with a spectral bandwidth set by the tunable filter (filter used here is 350 GHz). For these data Δf_(r)=1 kHz, a single interferogram is acquired in a laboratory time of 1 ms (only the central ˜8 μs is shown here). The time axis is also given in effective time, defined by the pulse-to-pulse time offset, ΔT, between the signal and LO pulses. The strong signal at zero time, which has been normalized to unity, corresponds to the overlap of the LO and unperturbed source pulse. The effect of the gas is to generate a free-induction decay signal that extends to longer times. For the single trace (upper, light), this FID tail is barely above the noise level. However, the coherently averaged interferogram (dark, offset for clarity-1000 averages shown) has greatly improved SNR;

FIG. 2B is the same data as in FIG. 2A with the vertical axis expanded 25 times. The averaged data clearly resolves the free induction decay (FID) signal. This FID of the molecules appears as a pulses due to the “rotational recurrences”, “rotational revivals”, or “commensurate echoes” as the rotating molecules rephrase;

FIG. 2C is a complex frequency-domain response generated by a Fourier transform of the averaged trace. Spectral response (magnitude only) is plotted here for a single trace (light) and 1000 averages (dark). The rf frequency axis corresponds to the lab time axis and represents the actual frequency of the digitized signal. The optical frequency axis corresponds to the effective time axis and shows the actual optical frequency relative to the 1560 nm cw reference laser discussed later.

FIG. 3A is an experimental setup for spectroscopy of a sample with two femtosecond erbium-doped fiber lasers phase-locked together with a well-defined difference in repetition rates and phase-coherent optical carriers. Pulses from the source laser (or frequency comb) are coupled into a free space section where a beamsplitter divides the source pulse into a signal pulse that passes through a cell containing H¹³C¹⁴N gas and a reference pulse that circumvents it. Because of the different path length travelled by the signal and reference pulses, there is a beneficial 2.6 ns time separation between them. The signal and reference pulses are then recombined with each other, and then with the LO pulse train from the LO frequency comb on a polarizing beam splitter (PBS). The combined pulses are coupled back into a single-mode fiber optic, filtered by a scanning bandpass optical filter and detected on a balanced detector (BD). The additional AOM and switch in the LO path are used to select the offset frequency of the LO such that the detected rf multi-heterodyne signal between the LO and source combs does not fall at baseband or Nyquist.

FIG. 3B is a detected voltage signal, with the peak normalized to one, corresponding to the overlap of the LO pulse with either the reference or signal pulse. Voltages are digitized synchronously with the LO pulses. The reference pulse arrives 2.6 ns sooner than the signal pulse and the separate arrival of the two is seen in the two separate peaks in the detected voltage. Note that the two peaks are easily distinguished from one another. In the terminology of conventional FTS, these two peaks correspond to two separate centerbursts in the interferogram. The time axis is given in both laboratory time and effective time for a Δf_(r)=3.14 kHz. The ˜320 microseconds duration in laboratory time corresponds to the time to acquire the N₁˜32,000 LO pulses, spaced by ˜10 nsec. The ˜10 nsec duration in effective time corresponds to the same N₁˜32,000 LO pulses with a point spacing set to equal the ΔT˜300 fsec time offset between the LO and source pulses. The data trace shown corresponds to actual data, averaged over 10⁴ individual traces using coherent real-time averaging.

FIG. 4A is a setup used to stabilize the two frequency combs such that they are both phase-coherent with sub-Hertz residual linewidths and with a very well-defined difference in repetition rates, Δf_(r). The locking setup is most easily illustrated in the frequency domain picture. A pair of teeth from each comb is stabilized to a pair of cw reference lasers. The phase-locks use identical rf offsets for both combs. Repetition rates are chosen here such that there are exactly N₁ LO comb teeth between the two points and N₂=N₁−1 for the source comb. As a consequence, in the time domain, the interferograms repeat precisely every N₁−1 LO pulses.

FIG. 4B is a setup used to phase-lock a single comb to the two cw reference lasers. The two cw lasers are phase-locked to a high finesse optical cavity using a Pound Dreyer Hall lock. The stabilized light from the 1560 nm and 1535 nm lasers is combined on a fiber and enters on the left hand side of the figure as shown. PLL-Phase lock loop, ISO-optical isolator, PZT-piezoelectric transducer.

FIG. 5A is Data processing for the signal from a 15 cm long cell filled with 25 Torr of H¹³C¹⁴N gas in which the Time domain interferogram signal after 10,000 averages. The trace is divided into reference and signal portions as shown. The reference corresponds to the cross-correlation of the LO pulses with the source pulse that circumvented the sample while the signal corresponds to the source pulses that passed through the sample. (See FIG. 3A).

FIG. 5B is a Fourier transform of reference, {tilde over (V)}_(R)(v_(k)), and signal, {tilde over (V)}_(S)(v_(k)), sub-traces (magnitude only, phase not shown). The signal shows the absorption of about 6 ro-vibratonal HCN lines.

FIG. 5C is a normalized response, {tilde over (H)}_(λ)(v_(k)), in magnitude (dark) and phase (light) calculated from the ratio of the complex, frequency-domain signal and reference of FIG. 6 b after shifting the frequency axis by the offset frequency, v₀. Dotted line shows the reference trace with the same offset. A spurious peak at 5.29 THz due to rf pick up is visible here but removed from the final data. When response traces are concatenated to cover the full comb spectral width, only the data within the filter FWHM are used.

FIG. 6A is the measured complex spectral response, {tilde over (H)}(v), of the ro-vibrational lines corresponding to the first overtone H—C stretch vibration in HCN over a 9 THz bandwidth at 220 MHz resolution. The “p-branch” and “r-branch” are clearly visible. For each isolated ro-vibrational line, the magnitude corresponds to the standard gas absorption profile (divided by two since absorption is typically given for intensity) and the phase corresponds to the delay of the light due to the change in the index of refraction. Three different HCN gas pressures are shown: 25-Torr in a 15-cm long cell (light), 2.7-Torr in a 20-cm long cell (dark), and an empty cell. Each trace is generated by concatenating the complex frequency response (e.g. FIG. 5C) measured at forty five different settings of the tunable 2 nm-wide optical bandpass filter.

FIG. 6B is an expanded view of the complex HCN spectrum. For the magnitude data, the three traces have been offset for clarity. For the higher-pressure cell (25 Torr), the magnitude has the Lorentzian profile expected from collisional line broadening with a corresponding “derivative”-like shape for the phase profile. For the lower-pressure cell (2.7 Torr), the magnitude has a gaussian profile expected from Doppler line broadening, again with corresponding “derivative”-like shape for the phase profile. The molecular hot bands can clearly be seen as much smaller magnitude and phase excursions. The combined benefits of extended averaging and normalization by the reference channel can be seen in the flatness of the baseline. The statistical noise averaged over the full spectrum is ˜0.04% in magnitude and 400 μrad in phase. Over the central portion of the spectrum, it drops to ˜0.025% in magnitude and 250 μrad in phase. These data were acquired with 60 seconds of averaging at each of the forty five positions of the tunable optical filter giving a total acquisition time of 2700 sec. Of course, the SNR scales with the square root of time so that a 270 sec total acquisition time reduces the SNR in magnitude and phase by a factor of only ˜3.3.

FIG. 7A is the reconstructed time-domain response and joint time-frequency domain response of HCN gas corresponding to the complex frequency response of FIG. 6. This different, but equivalent, view of the data allows for a different physical picture of the gas response as discussed in the text. FIG. 7A is the total time-domain response over a 6 ns window at 55 fs point spacing (9 THz bandwidth) with the peak at t=0 normalized to unity. The time-domain SNR is 0.55×10⁶.

FIG. 7B is the expanded view of the time-domain response for the 2.7-Torr cell (dark), 25-Torr cell (light) and empty cell (gray). While the empty cell trace is flat, the filled cells clearly show a tail from the FID of the excited HCN molecules. For the lower pressure, 2.7-Torr cell, 170 rotational recurrences are visible out beyond 1.8 ns. Because of causality, the FID appears only at positive times (where t=0 corresponds to the arrival of the unperturbed excitation pulse). This single-sided interferogram is a direct consequence of the fact that the LO pulse train passes around rather than through the HCN sample. The envelope of the overall FID signature decays more slowly for the lower pressure cell because of the much lower collisional decoherence rate.

FIG. 7C is a further expanded view where the 25-Torr and 2.7-Torr cell data are offset for clarity. The low noise level can be seen in the middle empty cell trace (gray). The complicated time-domain structure is a result of interferences between the different ro-vibrational levels.

FIG. 7D is a Sonogram (short-time Fourier transform) of the 2.7 Torr data processed with 19 GHz frequency resolution and 52 ps time resolution. The bright vertical stripe at t=0 corresponds to the arrival of the signal pulse and has been normalized to 0 dB. At this frequency resolution the trailing decay signal from each absorption line is clearly visible. At 2.7 Torr, decay is driven primarily by Doppler dephasing, allowing for relatively long decay times.

FIG. 7E is a Sonogram under identical conditions as FIG. 7D but for the 25 Torr data. At this pressure collisions accelerate the decay time and collisional resonance between like states is evident in the thermal distribution appearance to the decay times.

FIG. 8A is the Coherent dual-comb spectrometer reconfigured from FIG. 3A to instead perform spectroscopy on arbitrary cw waveforms, in other words to measure the spectrum of the cw waveform versus time. As before, the two combs are femtosecond Er-fiber lasers optically phase-locked together with repetition rates, f_(r)≈100 MHz differing by Δf_(r). A frequency-shifted duplicate of comb 2 is generated with an acousto-optic modulator (AOM). The comb sources output a set of coherent pulses that serve to down-sample the cw laser's E-field. The resulting voltages are digitized synchronously with comb 1, recorded and processed. For spectroscopy of a sample, we referred to the two combs as a source and LO comb. Here the combs both act on equal footing to sample the cw laser field and we refer to them as comb 1 and comb 2. The system is clocked synchronously with comb 1.

FIG. 8B is the wideband spectrum calculated from the product of both sampled voltages, as discussed below, to provide the absolute frequency calibration with a time resolution Δf_(r) ⁻¹=T_(W) (320 μs and 30 μs).

FIG. 8C is a High-resolution measurement of the instantaneous laser frequency (right) and corresponding high-resolution spectrum (left). Measurements at multiple contiguous or spaced time intervals T_(W) can be concatenated to monitor the laser frequency output over long times.

FIG. 9A illustrates the sampling of the cw laser E-field by pulses from the Combs 1 and 2. The pulses arrive at intervals of T₁=f_(r1) ⁻¹ and T₂=f_(r2) ⁻¹, that differ by ΔT=T₁−T₂. (Not to scale).

FIG. 9B is the real part of the sampled voltages V₁ and V₂ of the laser E-field by comb 1 and comb 2 for zero ceo phase.

FIG. 9C is the real part of the complex product V₁V₂* rescaled to a point spacing equaling ΔT to correspond to the incremental time delay between the sample pairs V₁ and V₂.

FIG. 9D is the Frequency-domain picture of the input cw laser field (grey), comb source 1 (dark) and 2 (light). Zero on the frequency axis corresponds to the teeth of combs 1 and 2 phase-locked to the underlying cw reference laser.

FIG. 9E is the instantaneous frequencies f_(1,2) calculated from V_(1,2), which equal the laser frequency with an ambiguity nf_(r1,2).

FIG. 9F is the wideband spectrum calculated from the Fourier transform of V₁V₂*, which resolves this ambiguity.

FIG. 10A is the mean frequency of a external cavity diode laser (ECDL) where the measurements are spaced by 1.6 ms. The mean frequency is calculated from the wideband spectrum. The steep slopes correspond to tuning of the ECDL grating at maximum speed, and are corrected for coherence-induced biases. The darker color indicates required a priori knowledge of the laser frequency of ΔT⁻¹=3.2 THz (25 nm).

FIG. 10B is the wideband laser spectrum as the grating is stopped. See FIG. 11E for individual spectra.

FIG. 10C is the high-resolution spectrum (false color, log scale) for the region indicated in FIG. B, normalized to unity area, with the absolute frequency independently assigned from the simultaneously acquired wideband spectrum. The darker blue indicates the Nyquist bandwidth of f_(r)/2.

FIG. 10D is the instantaneous laser frequency over T_(W) (top panel) and the corresponding high-resolution spectra (bottom panel), for three specific measurement times t₀. The frequency offset is 194.423000000 THz±1 kHz

FIG. 10E is an example wideband spectra for times when the laser frequency is increasing (at 1.2874 s and 1.3433 s), and when it is relatively quiet (at 1.5714 s). The spectra are normalized to a 0 dB peak. The measured signal is the convolution of the intensity spectrum and the effective instrument response function, which is determined by the residual phase noise between the combs. For the trace at 1.5714 s, the cw laser frequency fluctuations are well below 100 MHz and ideally one expects a delta function. The SNR for this trace is ˜15 dBc as measured from the peak to the close-in noise pedestal, or ˜50 dBc/Hz when normalized by the 3.13 kHz resolution. At larger frequency offsets, the SNR improves to 26 dBc, or 60 dBc/Hz. This SNR is given by the magnitude of the relative phase noise between the combs, as opposed to the magnitude squared. Therefore, low residual phase noise (and linewidth) between the combs is particularly critical in narrowing the instrument response function and achieving a reasonable SNR. In addition, there can be a broadband noise contribution from the unwanted “sum frequency signal”. For the data at 1.5714 s, this sum frequency beat appears offscale but is easily identified by its broader width and ignored. For the other time points, the sum frequency signal is roughly equally distributed over all Fourier frequencies and causes the raised noise floor. (IQ detection could remove this component in principle and bring these noise floors down to the level exhibited for the 1.5714 s data.).

FIG. 11A is the mean frequency of a laser mechanically perturbed at t≈2.5 ms continuously monitored at a time-resolution of T_(W)≈30 μs. When the laser frequency changes by more than 50 MHz over T_(W), the mean optical frequency is calculated from the wideband spectrum and corrected for coherence-induced biases. When the laser frequency changes more slowly, the mean optical frequency is calculated from the high-resolution (34 kHz) spectrum. The 20 GHz fluctuations are within the ΔT⁻¹≈0.15 THz ambiguity. The cw laser frequency is 192.293980842 THz±1 kHz at t=0, but could fall anywhere within the 5 THz comb bandwidth, unlike cw laser heterodyne spectroscopy where the source and reference lasers must closely match.

FIG. 11B is a power spectral density of FIG. 11B over a time window extending out to 0.18 s showing the mechanical system resonance at ˜1 kHz. The frequency noise floor is a result of the actual ˜10 MHz scale technical fluctuations in the laser frequency.

DETAILED DESCRIPTION

We first give an overview of coherent dual-comb spectroscopy of passive samples and then an overview of coherent dual-comb spectroscopy applied to active cw sources (beginning in paragraph 0093).

With reference to FIG. 1A, the basic approach of coherent dual comb spectroscopy is to interfere two combs with repetition periods differing slightly by Δf_(r). If the two combs are combined before the sample, the system is analogous to traditional FTS and yields the intensity absorption spectrum only. Alternatively, the source comb passes through the sample and afterwards is combined with a local oscillator (LO) comb. This approach is analogous to asymmetric or dispersive Fourier transform spectroscopy (DFTS) and yields both the magnitude and phase spectrum of the sample. Note that for the single spatial mode output of combs, alignment issues are potentially less troublesome than in DFTS.

The operation of the dual comb spectrometer can be viewed in either the time domain or frequency domain (FIGS. 1B, 1C and 1D). In the time domain, the pulses from each comb source overlap on the photo detector at varying time delays, essentially creating a virtual scanning interferometer as in asynchronous optical sampling. The resulting cross-correlation of the source and LO pulses is analogous to a conventional interferogram in DFTS.

An example interferogram is shown in FIGS. 2A-2C with a centerburst followed by a trailing electric field, which, in the case of a molecular gas sample considered here, is just the free-induction decay (FID) of the molecules. A Fourier transform of this interferogram yields the spectral response of the gas and, since the interferogram is single-sided, one obtains both magnitude and phase. Because the pulse trains are periodic and repeatedly move through each other, there is no delay between interferograms. Rather a new interferogram begins to form as soon as the previous is finished. It should be understood that the down sampling nature of this experiment introduces two time scales, which will be referred to herein as “lab time” and “effective time.”

Lab time is the time taken to construct an interferogram and is characterized by a point spacing equal to the comb repetition period. Effective time corresponds to the time scale of a single reconstructed pulse and is characterized by the much smaller time step equal to the difference in comb repetition periods, given by ΔT=Δf_(r)/f_(r) ², where f_(r) is the comb repetition rate (FIGS. 1B, 1C). Nyquist sampling limits require that the instantaneous optical bandwidth of the combs satisfy the relationship Δv_(comb)<(2ΔT)⁻¹. Larger optical bandwidths can be measured by sequential measurements at different center frequencies, as done here.

In the frequency domain (FIG. 1D), the system is essentially a massively parallel laser heterodyne spectrometer. The phase and magnitude of each tooth of the source comb is measured through its heterodyne beat against one tooth of the LO comb. In order to insure a one-to-one mapping of the rf comb to the optical teeth, we require Δv_(comb)/f_(r)<f_(r)/(2Δf_(r)), which is equivalent to the Nyquist sampling condition, Δv_(comb)<(2ΔT)⁻¹ arrived at in the time-domain picture. Based on this Nyquist sampling constraints, for a given comb repetition rate and optical bandwidth, the minimum required time to acquire a spectrum equals Δf_(r) ⁻¹.

With this simple picture, several important questions are immediately apparent: how to normalize the signal to remove the distortions from the LO and source spectra, how to achieve high signal-to-noise ratio (SNR) at data rates of f_(r)˜100 MHz or more (i.e. much higher than conventional FTS), and what are the critical requirements on the frequency combs.

Normalization: The source light that passes through the sample is only measured with respect to the LO comb light and both can contain significant spectral structure in addition to etalon effects occurring in the beam path. Moreover this structure can vary with time. Without a rapid normalization scheme to remove these distortions, the effective SNR of the trace will be completely limited by background structure. This may be accomplished with a second interferometer. Here, a portion of the available interferogram time window is used by time-multiplexing a signal and reference pulse (FIGS. 3A, 3B and 5) to allow an update to the reference every 300-1000 μs (Δf_(r) ⁻¹). This approach sacrifices the potential 100-MHz resolution set by the comb repetition rate for the ability to measure a fully normalized spectrum in a compact setup with a single detection channel, which is particularly well suited to the single-sided interferogram because of the larger available time window. Because the reference and signal pulses travel through nearly identical optics etaloning effects are also strongly suppressed, allowing for an extremely flat baseline. An additional advantage of normalizing the data on a sub-millisecond time scale is the suppression of 1/f noise on the laser sources.

Signal to Noise Ratio: The single point detector shown in FIG. 1 has the benefit of simplicity, high frequency accuracy, resolution, and adaptability to the mid IR and far IR. However it may suffer an intrinsic SNR limitation compared to other comb spectroscopy that uses array detection because of detector dynamic range, and broadband laser technical noise. The high frequency resolution and accuracy of a dual comb spectrometer are not of particular use unless the spectrum is measured with sufficient SNR. One route to higher absolute signal is a longer interaction length using a cavity or multi-pass cell. Cavity-enhanced single comb spectroscopy has been well demonstrated and cavity-enhanced dual comb spectroscopy has been recently demonstrated at 1 micron.

Regardless of the signal strength, one would desire as high a SNR as possible. The SNR can be improved by averaging over many interferograms in analogy with co-adding of spectra from rapid-scanning FTS systems. However, given the high digitizer rates (set by the repetition rates of ˜100 MHz), the acquisition of multiple interferograms will lead to very cumbersome data file sizes, and furthermore individual phase-correction of these interferograms would require unacceptably long processing times. A solution for both these problems is to implement real-time coherent averaging by forcing the interferograms to have identical phase and summing them in real time. This coherent averaging can easily reduce file sizes by a factor of 10⁴ or more. Coherent averaging is possible only while the two combs remain phase coherent. Roughly speaking, coherent averaging requires a relative linewidth between the combs equal to the inverse of the averaging time, i.e. a relative linewidth much narrower than Δf_(r).

Frequency Comb Stabilization: FIGS. 1A-1E assume ideal source and LO frequency combs, however it is the real-world spectral and noise properties of the combs and their associated stabilization schemes that determine the frequency accuracy, resolution (and spectral coverage) of the spectrometer.

A frequency comb by itself has no inherent accuracy, as the free-running comb's two degrees of freedom can drift randomly; frequency accuracy is only obtained if the comb is either locked to or monitored against some pair of absolute references. For example, if one tooth of the comb is locked to an underlying cw reference laser, as is done here, the frequency of each tooth can then be inferred from this cw reference laser and the counted repetition rate. Without an underlying reference, accuracy must be derived from a known absorption line. The resulting accuracy is typically several orders of magnitude worse than is possible with stabilized combs.

The frequency resolution is in principle limited by the comb repetition rate. However, a resolution equal to the repetition rate can only be achieved if there is no apodization and if the residual phase noise between the combs is low.

As seen in FIG. 1E, a resolution equal to the comb repetition rate requires a residual comb linewidth of Δf_(r). (The residual linewidth is defined as the measured linewidth between teeth of the two combs, the absolute linewidths can be much larger). Otherwise, comb teeth in FIG. 1E will overlap ruining the one-to-one correspondence between the optical comb lines and rf comb lines, and limiting resolution to >f_(r).

Residual linewidths below Δf_(r) can enable coherent averaging, as described above, and can also allow a resolution below f_(r). Strictly speaking, the system only samples at the discrete comb tooth frequencies, but one can in principle shift the combs and repeat the measurements to acquire data at interleaved comb frequencies. However, since spectral linewidths narrower than the repetition rate will have a time domain structure that extends across multiple interferograms, achieving resolution below f_(r) requires continuous (unapodized) acquisition and processing of multiple contiguous interferograms. This, in turn, requires mutual coherence between the combs across the multiple interferograms. In other words, to have a resolution n times better than f_(r) one must have residual linewidths below Δf_(r)/n

Within the context of the above discussion, multiple approaches to stabilizing the frequency combs have been used in dual comb spectroscopy. The use of free-running combs, i.e. femtosecond lasers, is certainly the simplest method and was used in a pioneering demonstration of dual comb spectroscopy; however such an approach provides no absolute frequency accuracy, has limited frequency resolution and coherent averaging is challenging (although coherent averaging has been achieved in the far-IR where carrier phase is less of a problem). In an intermediate approach, the carrier frequency of each comb is phase-locked to a common reference cw laser, and the pulse repetition rates are monitored. The carrier frequency of the free-running combs may alternatively be monitored and the signal corrected in post processing. This approach removes the main source of drift—the carrier phase—and can provide high frequency accuracy if the reference laser frequency is known. To further enhance the coherence between the combs both combs can be fully phase-locked to a pair of calibrated cw reference lasers (the approach taken here). The accuracy of the combs is set by the cw reference laser frequency, which is measured with a separate self-referenced comb against a Hydrogen maser and is accurate to ˜10 kHz (limited by drift in the cw reference between measurements). This accuracy generally far exceeds the statistical uncertainty in the line center of a Doppler-broadened line (even for the high SNR signals demonstrated here) and a looser cw reference with ˜MHz accuracy locked to a simpler cavity, molecular line or even fiber loop would suffice. With our approach, the phase locks are sufficiently tight such that the residual linewidths are ˜0.3 Hz and allow for long coherent averaging times of 3 seconds and a corresponding improvement in SNR.

If high phase stability is required over a very broad bandwidth, one could self-reference the comb using the standard carrier envelope stabilization techniques in combination with a lock to a single optical reference.

The optical layout of the experimental system is described next. With reference to FIG. 3A, the two comb sources are erbium-doped femtosecond fiber lasers with repetition rates f_(r,LO)≈f_(r,S)≈100 MHz that differ by Δf_(r)=3.14 kHz for most of the data shown here and are centered around 1560 nm. Δf_(r) yields the effective time step ΔT=Δf_(r)/(f_(r,S)f_(r,LO))between the source and LO pulses. The comb outputs are amplified to 20 mW in a short erbium-doped fiber amplifier which also provides a small amount of spectral broadening.

A splitter after the amplifier (not shown in FIG. 3A) directs a portion of the light to the phase-locking setup described below. The remainder of the source pulse train is directed through the interferometer signal and reference paths and polarization multiplexed with the LO pulse train. The combined reference and signal then pass through a tunable grating monochromator which filters the light prior to detection and digitization. For data shown in FIG. 6-7 the filter has 2 nm (200 GHz) full-width half-maximum (FWHM) Gaussian profile. Data taken with different filter center wavelengths are separately processed, as described later, and then coherently stitched together to generate the total spectrum.

The optical filter is not necessary and an alternative design is to use the full comb spectrum. On one hand, the filter reduces the power on the detector and requires multiple measurements to reconstruct the spectrum. However, in a shot noise limited system this additional measurement time is exactly offset by the lower noise of a reduced spectrum. Moreover, in the dynamic-range limit (set by detectors or digitizers), filtering the spectrum can lead to an overall increase in SNR. Frequently, the system is dynamic range limited near the center of the spectrum and the filter improves the SNR.

Spectral filtering allows the system to simultaneously meet Nyquist sampling conditions and suppress 1/f noise. Specifically, the filtered comb bandwidth (Δv_(comb)) must be limited to Δv_(comb)≦1/(2ΔT)=f_(r,S)f_(r,LO)/(2Δf_(r)), to avoid aliasing. Thus a broad spectral bandwidth requires a low value for the difference in comb repetition rates, Δf_(r). However, the interferogram updates at Δf_(r) and will be degraded by 1/f noise in the comb spectra at low values of Δf_(r). The use of the filter allows us to maintain the high scan speed (Δf_(r)) by limiting the instantaneous spectral bandwidth, and later recover the entire spectrum through coherent stitching.

In the disclosed non-limiting embodiment, after the filter, the orthogonally polarized signal and LO pulses are mixed with a fiber-coupled polarizing beam combiner oriented at 45° with respect to the laser polarization states. Polarization paddles are used to fine tune the balancing of light levels. The two outputs of the beam combiner are incident on a 110 MHz commercial balanced detector. Balanced detection suppresses the strong homodyne signal from the individual pulses and allows maximization of the dynamic range available for the heterodyning between the pulses. The output of the detector is low pass filtered at 50 MHz to avoid aliasing of signals above the Nyquist frequency and to relax the timing between the digitizer clock and the pulse arrival on the photodiode. The detector signal is then digitized by a 100 MHz 12-bit digitizer.

As shown in FIG. 3A, the LO passes through a 120 MHz acousto-optic modulator (AOM) and fiber switch to select either the first (shifted) or zeroth (unshifted) order beam from the AOM. There are two benefits to the AOM. First the LO power can be adjusted via the AOM's rf power to avoid photodetector saturation as the optical filter is tuned from high to low power regions of the comb spectra. Secondly, at certain optical filter positions the rf beat between the combs falls near the Nyquist frequency (f_(r,LO)/2) or near 0 Hz. At either of these frequencies, aliasing effects lead to spurious data. These regions are avoided by shifting the rf interferogram using the first-order deflected beam thus shifting the LO by ˜f_(r,S)+f_(r,S)/4 (˜125 MHz). Maintaining carrier phase of the interferogram during coherent averaging does require the AOM be driven at an integer multiple of Δf_(r), but this shifting is otherwise straightforward. One could suppress these reflections about Nyquist with IQ detection however it is difficult to realize more than 20 dB of suppression from an optical IQ detector and 1/∫ noise would remain near 0 Hz.

With reference to FIGS. 4A-4B, the basic method to phase-lock the two Er fiber combs together are schematically illustrated. A pair of comb teeth from each comb are locked to two cw reference lasers at 1535 nm and 1560 nm respectively. With this setup, in the frequency domain the source and LO combs exactly overlap every K and K+1 modes respectively, where K+1=N₁, defined earlier. In the time domain, the LO and source pulses exactly overlap every K pulses and K+1 pulses respectively. K sets both the number of samples in a single interferogram and the scaling between lab time and effective time (e.g. FIG. 3B axes). For most data shown here, K=31834. The integer length of the interferogram greatly simplifies coherent signal averaging as discussed below.

Both cw reference lasers are stabilized to the same high-finesse optical cavity via a Pound Dreyer Hall lock (PDH). In practice, the measured frequency of the 1560-nm laser and the measured repetition rates of the two combs, with respect to a Hydrogen maser are used to calculate the optical frequency of the source and LO comb teeth. Accuracy is linked directly back to the Hydrogen maser in this fashion. Note that in this case, the 1535 nm laser does not need to be cavity-stabilized as well but it is convenient.

In the disclosed non-limiting embodiment, the two combs are phase-locked to the cw lasers in an identical fashion. The comb light is combined with the two cavity-stabilized cw lasers, spectrally filtered to two channels at 1535 nm and 1560 nm with a 1 nm width, and detected on two 100 MHz detectors, as shown in FIG. 4B. Each detector detects a beat signal between the comb and one cw laser. The two beat signals are then used in two separate servos to tightly phase-lock the combs. Coarsely speaking, the 1560 nm light is used to stabilize the comb carrier frequencies and the 1535 nm light to stabilize the comb repetition rates (in truth these parameters are coupled). Because the optical bandwidth is relatively small compared to the carrier frequency, the carrier stabilization lock is the more critical of the two locks. For this lock, the rf beat signal is divided between the comb and 1560-nm laser by eight to increase the capture range and lock the divided signal to a 10 MHz reference by feeding back to an AOM with ˜100 kHz bandwidth and a much slower, but higher dynamic range, intracavity PZT . The servo to the PZT effectively removes much of the lower frequency noise on the repetition rate of the comb as well. However, to fully stabilize the repetition rate, we use the beat signal between the comb and the 1535 nm light, which is fed into a servo filter whose output is fed to the pump current for the laser cavity. By controlling the pump current, the repetition rate is controlled while roughly holding the carrier frequency fixed. The feedback bandwidth is ˜10 kHz limited by the laser dynamics combined with the erbium response. The 10 kHz is sufficient, however, phase lead compensation could achieve higher bandwidths.

The data acquisition is described next. The signal is digitized at the LO repetition rate to match the effective optical sampling of the LO pulses, which translates to a high sample rate. To reach a high SNR requires signal-averaging of multiple interferograms just as in conventional FTS. However, the brute force approach of directly digitizing the signal for as long as possible leads to massive data sets and makes deep averaging impractical. Therefore, instead coherent averaging is implemented by phase-locking the two combs as described above such that each interferogram is exactly an integer number (K) of points with exactly the same phase. In that case, successive interferograms in real time may be summed either onboard an FPGA or in software for a time approximately equal to the inverse of the relative linewidth between the combs. The requirement on the relative comb linewidth for coherent averaging is very demanding; we require linewidths <Δf_(r)/N_(A), where N_(A) is the number of coherently averaged interferograms.

In the disclosed non-limiting embodiment, the system coherently averages for ˜3 seconds corresponding to N_(A)=10,000 interferograms and a ˜0.3 Hz relative comb linewidth. The commercial digitizer must be triggered for each successive interferogram. A trigger, accurate to below the 10 ns clock period, is achieved with a digital pulse counter to divide the LO clock signal by K (FIG. 3A). The digitizer may require a small (˜20 μs) dead time between triggers which leads to a slight loss in frequency resolution (˜7%). At a 200 Mbyte per second sampling rate, a 600 MB data stream is reduced to 60 kB (or by 10⁴) with a corresponding decrease in processing time because the FFT is performed on the shorter data set. For acquisition periods greater than 3 seconds, drifts in optical path lengths and noise in the phase locks leads to a relative carrier phase drift between combs. Therefore, the sequential 3-second coherently summed interferograms must be phase corrected, based on the phase of the centerburst corresponding to the reference peak, and then sum them to generate an interferogram at longer averaging times.

Phase correction may also be applied more frequently to relax the requirements on the mutual coherence (residual linewidth) between combs but the processing load and data transfer from the digitizer will add additional dead time. However, more sophisticated FPGA's could make this realtime correction practical.

The first step of the data processing is to change the scaling of the data from laboratory time to effective time. Given the locking scheme used here, this scaling is straightforward. The sample point spacing is set to equal the effective time step ΔT=Δf_(r)/(f_(r,LO)f_(r,S))=1/(Kf_(r,S)). The three quantities Δf_(r), f_(r,LO), f_(r,S), are measured with frequency counters.

The time multiplexing of reference and signal pulses discussed earlier yields the optical interferogram, V(t_(i)) seen in FIG. 5A, where t_(i) is in effective time, with the two pulses separated by 2.6 ns and the signal peak is defined to be zero time. A “reference-only” time trace is generated by gating out all but the reference pulse, to form V_(R)(t_(i))=<V(−4.7 ns<t_(i)<−1.3 nsec)> and zero-padded to a length K. Similarly, a “signal-only” time trace is generated by gating out all but the signal pulse to form V_(S)(t_(i))=<V(−1.3 nsec<t_(i)<4.6 nsec)>, which is also zero-padded to K points. To deconvolve the reference from the signal, the two signals are first Fourier transformed to yield their respective complex spectra (FIG. 5B). The ratio yields the complex molecular response, shown in FIG. 5C, for the filter centered at λ as {tilde over (H)}_(λ)(v_(k)))={tilde over (V)}_(S)(v_(k))/{tilde over (V)}_(R)(v_(k)), where the tilde represent a Fourier transform evaluated on the optical frequency grid v_(k) with a known optical offset v₀ discussed below.

Because the time trace is padded to K points, the frequencies v_(k) correspond precisely to comb teeth of the signal comb. The resolution of the system is set by the temporal length (i.e. apodization) of V_(S)(t) giving a resolution of 1/4.6 ns≈220 MHz. The reference pulse is more strongly apodized with less spectral resolution to limit the contribution of noise from the reference to the final deconvoluted signal. Generally the reference pulse window must only be broad enough to contain any etaloning features (<1 ns) that also appear on the signal.

The optical frequency offset, v₀, is calculated from the measured frequency of the 1560 nm cw laser, the known rf offsets in the phaselocking, the known AOM frequency shift if present, and the measured repetition rate. For example, for the unshifted LO comb, v₀=v₁₅₆₀+f_(lock)+pKf_(r,S) where v₁₅₆₀ is the measured frequency of the 1560-nm cw reference laser, f_(lock), is the rf offset of the phase lock to the 1560 nm laser, and p is an integer chosen such that v₀ is close to the optical bandpass filter center. The term pKf_(r,S) reflects the ˜1.5 THz Nyquist ambiguity resulting from the ΔT time step. In addition the sign of the frequency axis of spectra in FIG. 5B is changed from positive to negative as necessary. Assigning the negative sign and values for p can be done with ˜1 nm level knowledge of the filter position. Thus applying these shifts correctly requires only a coarsely calibrated optical filter.

Despite the normalization, some slow background wander may remain on {tilde over (H)}_(λ)(v_(k)) due to multiplicative phase noise. While small (˜0.1% in magnitude or 1 mrad in phase), this wander can lead to ripple on the concatenated data. Therefore, we first remove it by separately fitting the magnitude and phase profiles of {tilde over (H)}_(λ)(v_(k)) to a third-order polynomial over a range covering the FWHM of the optical filter. The magnitude is then normalized to unity by dividing by the magnitude fit, and the average phase is set to zero by subtracting the phase fit. Strong absorption features can throw off the fits, so we mask out features that deviate from the baseline by more than 2%. For the phase profile it is sufficient to apply a 7 GHz box-smooth prior to fitting. Note that this filtering technique would have to be reconsidered if one were looking at very broadband spectral features.

The individual {tilde over (H)}_(λ)(v_(k)) are then coherently stitched together, i.e., concatenated, to generate an overall response function covering the full frequency span and at a spacing f_(r,S). Because the combs are phase-locked to stabilized cw lasers, the optical frequencies are stable over the measurement period and there are no issues in the concatenation. This processing involves only fast Fourier transforms, which are indeed fast if K is factorable into small prime numbers.

The measured sample response can be viewed in the frequency, time, and joint time-frequency domains. With reference to FIGS. 6A and 6B, the measured frequency-domain response is in terms of both transmission and phase shifts for the ro-vibrational band corresponding to the first overtone C—H stretch in hydrogen cyanide (HCN) for a 2.7-Torr cell and a 25-Torr cell, along with an empty reference cell. The frequency domain response is given by {tilde over (H)}(v)=1+4 π²ic⁻¹vχ(v)L, where L is the length, χ the linear susceptibility of the sample, c is the speed of light, and v is the optical frequency. For weak gas absorption {tilde over (H)}(v)≈1−α(v)L/2+iΔk(v)L, where α(v) is the usual absorption coefficient and Δk(v) is the corresponding phase shift. For the case of simple, well separated absorption lines considered here, the absorption lines appear with the standard Voigt profile and the phase shift with the corresponding “derivative”-shape. (These are actually just the real and imaginary parts of the complex Voigt function or error function for complex arguments.) Note that this signal differs from the simple absorption spectrum, α(v), usually generated in grating spectrometers or typical double-sided FTS, because we retrieve both the magnitude, α(v)/2, and the phase. Of course, using Kramers-Kronig relations one can in principle retrieve the phase from the magnitude spectrum; by measuring it directly any complications due to finite acquisition bandwidth are avoided.

In the time domain (see FIGS. 7A-7C), the signal is simply the Fourier transform of the above frequency-domain signal convolved with a sinc function (corresponding to the spectral bandwidth, Δv, of the normalized response) or H(t)≈e^(−i2πv) ^(c) ^(t)sinc(πtΔv)

[δ(t)+4π²ic⁻¹v_(c){tilde over (χ)}(t)L], where the Fourier transform of vχ(v) is approximated as v_(c){tilde over (χ)}(t) where v_(c) is the average carrier frequency. Physically, the forward scattered light or free induction decay of the excited molecules are observed. The initial pulse excites a superposition of ro-vibrational states. These vibrating states act as dipole emitters. However, they quickly rotate out of phase with each other to give no coherent forward scattered signal. Because the rotation rates are quantized, they rotate back in phase a short time later (set by the rotational constant), giving rise to another burst of coherent radiation in the forward direction, long after the excitation pulse has ended (FIGS. 2B or 7C). The mathematical description follows earlier work on two-level systems by Brewer and coworkers, or the work on commensurate echoes in THz time-domain spectroscopy or rotational recurrences in pump probe spectroscopy.

Of course, there is no particular reason to view the signature in the purely time or frequency domains, and in FIGS. 7D-7E a joint time-frequency domain picture is shown at two different gas pressures. In this picture, the band at zero time is the reconstructed initial pulse and the trailing signal at later times is the frequency-resolved free induction decay amplitude. Each HCN absorption line is (slightly counter-intuitively) evident as a trailing signal. The effect is due to the slowing of the light near the resonance which delays the transmitted light at that frequency. The overall signal decays as a result of collisions and Doppler effects. The collisions cause a phase-shift or de-excitation of the vibrating, rotating molecules so that they no longer add coherently. The Doppler shifts can be viewed as the molecules moving out of their original position so that their emitted dipole radiation no longer adds phase-coherently in the forward direction. Collisional effects dominate for the higher pressure cell and the effects of resonant collisions is evident in the faster relaxation of the more highly populated ro-vibrational states. Doppler relaxation dominates for the lower pressure cell.

The nominal signal-to-noise ratio will depend on whether it is quoted in the time or frequency (or joint time-frequency domains). The time-domain SNR is measured as the ratio of the peak of the interferogram to the standard deviation of the noise as measured from either the empty cell or, equivalently, before the centerburst. At a point spacing of 55 fs over the 6 ns time window, we obtain a time-domain SNR of

SNR_(t)=1.1×10⁴T^(1/2)   (1)

as a function of total acquisition time, T. For the longest acquisition time of T=45 minutes for the data of FIG. 8, we measure SNR_(t)=0.55×10⁶. (Note that in calculating this SNR we use the noise near the centerburst, which is larger by √{square root over (2)} than the noise at longer time offsets due to the effects of deconvolution.) This noise carries over to the frequency domain, where we measure an average SNR across the 41,000 frequency elements that span the full 9 THz of

SNR_(f)=50 T^(1/2)   (2)

in radians or fractional magnitude change in the reconstructed spectrum. The peak SNR near the center of the spectrum is about twice as high or 100 T^(1/2) over a 1 THz window.

If the time-domain noise is white, then the transformation from Eq. (1) to Eq. (2) is straightforward. Specifically, SNR_(f)=SNR_(t)/√{square root over (M)}, where M=9 THz/220 MHz=41,000 is the number of resolved frequency points. Although the frequency noise is not completely white due to variations in the comb power across the spectrum, this relationship roughly holds because SNR_(t)/√{square root over (M)}=54 T^(−1/2) in agreement with Eq. (2).

In terms of sensitivity for trace gas detection, the minimum detectable absorption Lα_(min) is often used as system metric. Eq. (2) gives the SNR for the magnitude, Lα(v)/2, so that one might consider twice that value as a measure of the sensitivity. However, this relation oversimplifies the situation for several reasons. First, the measurement is made at much higher resolution than a typical spectral linewidth; the signal across the entire line contributes to the overall sensitivity and there is a corresponding increase in sensitivity over Eq. (2) of a factor of ˜√{square root over (Δv_(L)/v_(res))}, where Δv_(L) is the FWHM collisionally broadened linewidth and v_(res) is the system resolution. (This is the SNR enhancement for smoothing the data to the resolution of the spectra of interest). Moreover, the measurement is made over multiple spectral lines and the signal across all the spectral lines contributes to the overall sensitivity. Therefore a matched filter yields a sensitivity that is actually improved by a factor of √{square root over ((π/4)Σ_(i)α_(i) ²Δv_(Li)/(α₀ ²v_(res)))}{square root over ((π/4)Σ_(i)α_(i) ²Δv_(Li)/(α₀ ²v_(res)))} where the i^(th) line has peak absorption α_(i) and Lorentzain FWHM linewidth Δv_(Li) and the sum is over all ro-vibrational lines. With these extra factors, for HCN, the absorption sensitivity is improved over Eq. (2) by ˜10. However, even including these extra factors, the dual comb spectroscopy is hard pressed to compete in terms of sensitivity with a single swept cw laser spectrometer under identical conditions. The advantage of the dual comb system is the absolute frequency accuracy, high frequency resolution, and the broadband coverage that should allow for good discrimination against other absorbing species.

Coherent dual comb spectroscopy can be viewed as massively parallel heterodyne laser spectroscopy, infrared time-domain spectroscopy, or a form of dispersive Fourier transform spectroscopy. We provide a general discussion of some of the important attributes of dual comb spectroscopy including frequency accuracy, resolution, normalization, and sensitivity. These considerations motivated our particular dual comb spectroscopy design, which is described in detail. We demonstrate that reasonably high SNRs can be achieved with a single point detector over a large number of resolved frequency elements through coherent averaging. Detector arrays or selected filtering bring the possibility of more rapid acquisition of high SNR spectra. We note many of the SNR and accuracy considerations here apply as well to dual comb LIDAR. for further understanding of other aspects of the LIDAR and associated dual comb usage, attention is directed to U.S. patent application Ser. No. ______, entitled ______ which is assigned to the assignee of the instant disclosure and is hereby incorporated herein in its entirety. As a demonstration of the capabilities of dual comb spectroscopy, we measured the fully normalized, complex response of the ro-vibrational band of HCN for the first overtone of the C—H stretch. This measurement covers a 9 THz bandwidth with 41,000 resolution elements. The peak frequency domain uncertainty is 0.025% in magnitude giving a peak fractional absorption sensitivity of 0.05% (converting to intensity) and 250 microradians in phase.

There are many other potential configurations of dual comb spectrometers including the incorporation of remote fiber delivery, long interrogation paths, multipass cells or cavities, detector arrays etc. This approach may additionally be combined with multi-dimensional spectroscopy.

Another such configuration is a coherent dual-comb-based spectrometer reconfigured to measure an active cw source. Because the dual comb spectrometer can provide measurement with both high accuracy and high update rates, the system is capable of measuring cw optical waveforms at 30-μs and 320-μs time resolution over THz bandwidths. Within each time interval, the spectrometer returns the laser frequency spectrum with kHz absolute accuracy and time-bandwidth limited precision. Unlike etalon-based techniques, each measurement is independently calibrated, which allows for discontinuous source tuning between measurements and the characterization of arbitrary cw waveforms. Unlike previous dual-comb frequency metrology, this system is coherent so that the measurements can be made rapidly with time-bandwith limited resolution and accuracy. This configuration readily enables, for example, optimized waveforms for sensing applications including multi-species gas detection, coherent laser radar, and optical metrology.

Nearly all high sensitivity optical-based sensing relies on cw lasers because of their high photon flux, long interaction paths and ever-increasing frequency agility. Current cw lasers are capable of tuning rapidly or even hopping to different optical frequencies over THz of optical bandwidths, creating broadband cw optical waveforms, particularly in the telecommunication band. The challenge is that, unlike agile rf sources, the waveforms from free-running lasers are loosely controlled and not highly reproducible, particularly in optical frequency. Phase-locked laser systems have demonstrated well-controlled scanned or stepped optical waveforms using frequency combs, but with limited speed and/or bandwidth. For fast, wide-bandwidth measurements, a continuously tunable laser is usually swept in a quasi-linear ramp and referenced (i.e. linearized) against an etalon or, for higher linearity, phase locked to an interferometer. The dispersion and pathlength drift of the etalon/interferometer will limit the accuracy of these methods; these limitations can be overcome by referencing the swept laser to the evenly spaced modes of a self-referenced frequency comb for MHz accuracy. However, one is still restricted to a quasi-linear sweep with limited speed and any absolute frequency determination requires a separate measurement. Furthermore, the requirements on the cw laser source are restrictive as it must tune continuously since mode-hops or frequency jumps will destroy the calibration. More importantly, a linear sweep can be wasteful in acquisition time. For example in spectroscopy or optical metrology, the laser would ideally jump to spectroscopic features of interest and then scan slowly across the lineshape. Similarly, in LIDAR more complicated waveforms would enable optimized range-Doppler ambiguity functions or synthetic aperture LIDAR similar to RADAR systems.

A dual-comb spectrometer capable of the rapid absolute measurements of a dynamic cw laser thereby enabling the calibration of discontinuous/arbitrary cw waveforms from free-running lasers. Frequency combs, based on stabilized femtosecond lasers, have allowed for metrology of static lasers with unrivaled frequency accuracy.

It should be understood that a comb is an ideal ruler against which to measure a cw laser's instantaneous frequency, provided two issues are addressed. First, the phase noise contribution of the frequency comb must be negligible, which is achieved through sub-radian phase-locking to a stabilized cw reference laser. Second, the measurement ambiguity must be removed equal to the comb repetition rate. This can be done by a Vernier approach using combs with different repetition rates as suggested by Ma; Peng et al. have measured absolute laser frequencies within seconds in this way. Here this technique is extended to optically coherent combs, allowing for time-bandwidth limited measurements at sub-millisecond time resolution. The approach disclosed herein is similar to coherent dual-comb spectroscopy but applied to an active source rather than a passive sample.

With reference to FIG. 8A, an alternate non-limiting embodiment of the system is adapted from the coherent dual-comb spectrometer discussed above. The basic operation is outlined in FIGS. 9A-9F. The two comb sources have repetition rates f_(r1)≈f_(r2)≈f_(r)=100 MHz that differ by Δf_(r)=f_(r2)−f_(r1), and corresponding periods T₁=f_(r1) ⁻¹ and T₂=f_(r2) ⁻¹ that differ by ΔT=T₁−T₂. A comb source (e.g. comb 1) outputs a series of short pulses, E₁(t−mT₁)exp[iδφ₁(mT₁)] where m is the pulse number and for simplicity we assume zero carrier-envelope offset (ceo) phase so that δφ₁ represents carrier phase noise. These pulses are combined with the cw laser field, E_(L)(t)≡E_(L) exp[iθ_(L)(t)], with amplitude E_(L) and phase θ_(L)(t). We wish to determine the instantaneous frequency f_(L)(t)={dot over (θ)}_(L)(t)/(2π). Since the pulses are short, they act as a delta-function and their detected overlap with the laser field yields the sample voltages V₁(t₁)∝E_(L) exp(i[θ_(L)(t₁)−δφ₁(t₁)]) at discrete times t₁=mT₁. According to the sampling theorem, for an instantaneous laser linewidth below f_(r1)=T₁ ⁻¹, the instantaneous frequency of V₁(t₁) is f₁(t₁)=f_(L)(t₁)−δ{dot over (φ)}₁(t₁)/(2π)−nf_(r1), where the integer n reflects the ambiguity inherent in sampling non-baseband signals. In a comb picture, f₁ is the heterodyne frequency of the cw laser with the nearest comb tooth, labeled n. The above discussion equally applies to comb 2, resulting in V₂(t₂) at times t₂=mT₂ with instantaneous frequency f₂(t₂)=f_(L)(t₂)−δ{dot over (φ)}₂(t₂)/(2π)−nf_(r2). For negligible comb frequency noise, δ{dot over (φ)}_(1,2)→0, the measured f_(1,2) equals the instantaneous laser frequency f_(L), with an ambiguity nf_(r1,2), and the Fourier transform of V_(1,2) is a high-resolution laser spectrum with the same ambiguity.

The unresolved ambiguity nf_(r1,2) exceeds any a priori information on f_(L)(t), particularly for a rapid, discontinuously tuned laser. In the spirit of Refs. and Fourier spectroscopy to resolve this ambiguity, we consider the product V₁V₂*, which is a measure of the correlation of the E-field at differential sample times Δt=t₁−t₂=mΔT. Over a laboratory measurement time T_(W)≡Δf_(r) ⁻¹, this time coordinate Δt ranges from 0 to T₂ as the pulses of comb 1 raster through the pulses of comb 2. For negligible residual comb noise, i.e. (δφ₁−δφ₂)→0, and f_(L)(t) stable over T_(W), V₁V₂*∝|E_(L)|² exp(i2πf_(L)Δt) expressed as a function of Δt, and its Fourier transform with respect to Δt yields the “wideband” laser spectrum at resolution f_(r2)=T₂ ⁻¹ with a much coarser ambiguity ΔT⁻¹. In other words, we can now resolve the fine ambiguity, nf_(r1,2), to identify f_(L) to within ΔT⁻¹>>f_(r), which completes the measurement assuming some coarse a priori knowledge of f_(L). The total measurement bandwidth is limited only by the comb spectral width. We identify the required observation time of T_(W)=Δf_(r) ⁻¹ as our spectrometer's time resolution.

The above discussion assumed f_(L) changed slower than f_(r)/T_(W). However, this processing gives the wideband spectrum for broadband sources as well provided there is no coherence between spectral elements. A strongly chirped cw laser ({dot over (f)}_(L)>f_(r)/T_(W)) will exhibit such coherence and the spectrum is distorted, which must be corrected for in the frequency estimation. It should be understood that the operation of the spectrometer can be described in the frequency domain or time domain (see FIGS. 9D-9F and 9A-9C), where it is similar to a linear optical sampling system or Fourier spectroscopy

Complex signals V_(1,2), require IQ detection, but heterodyne detection yields only real signals. Consequently, the measured frequencies “fold” about a Nyquist frequency of f_(r)/2 for the high-resolution spectrum and (2ΔT)⁻¹ for the wideband spectrum. In particular, the wideband spectrum, which is calculated by multiplying V₁ and V₂, setting the point spacing to ΔT, and taking the Fourier transform is shown below to be proportional to the product I_(comb)(f_(opt)±nf_(r2)+qΔT⁻¹)I_(L)(f_(opt)±nf_(r2)+qΔT⁻¹), where I_(comb) is the comb intensity spectrum, I_(L) is the laser intensity spectrum, f_(opt) is an offset frequency equal to the frequency of the common comb tooth phase-locked to the underlying 1560 nm cw reference (or 1535 nm reference) laser, n runs from 0 to N/2, and q is an integer. In other words, we obtain the product of the comb and source intensity spectra evaluated on a frequency grid with spacing f_(r2), an offset frequency f_(opt), and Nyquist frequency (2ΔT)⁻¹. The spectral dependence of the comb can be removed by normalizing the signal with a simultaneous measurement of the total cw laser power. For a broad coherent source, the measured I_(L) will have distortions since the “measured” chirp is twice the actual chirp and this is taken into account in a frequency estimation. (See supplementary information).

Clearly, with heterodyne detection there are aliasing issues (i.e. reflections at the Nyquist frequency) for spectra near ∫=∫_(opt)+q(2ΔT)⁻¹, which is circumvented by also calculating the Fourier transform of the product of V₁ and V_(2S), the samples from the shifted comb 2, which is shown below to yield I_(comb)(f_(opt)+(4ΔT)⁻¹±nf_(r2)+qΔT⁻¹)I_(L)(f_(opt)+(4ΔT)⁻¹±nf_(r2)+qΔT⁻¹) assuming a shift 1.25 f_(r2); together the two measured wideband spectra remove the aliasing problems and allow assignment of the absolute frequency, on each measurement, to within (ΔT)⁻¹, regardless of the central carrier frequency. Similarly, the single heterodyne signals are combined to expand the ambiguity of the high-resolution measurement to f_(r2) rather than f_(r2)/2.

The choice of parameters is flexible. We use f_(r1)≈f_(r2)≈f_(r)=100 MHz and Δf_(r)=3.1 kHz and 34 kHz for FIGS. 10 and 11A respectively, with T_(W)=Δf_(r) ⁻¹≈320 μs and 30 μs and (2ΔT)⁻¹=f_(r1)f_(r2)/(2Δf_(r))≈1.6 THz and 0.15 THz, respectively. Note the tradeoff from T_(W)≡Δf_(r) ⁻¹≈f_(r) ⁻²ΔT⁻¹; finer time resolution, T_(W), implies a more restrictive ambiguity, ΔT⁻¹, and therefore requires better a priori knowledge of the laser frequency. The comb bandwidths are ˜50 nm to match the cw laser tuning range. The processing needs four Fourier transforms of f_(r)/Δf_(r)˜2¹⁵-point records, which requires ˜10 ms on a pc and is done in post-processing. However, field programmable gate arrays could provide real time output.

The two passively mode-locked fiber lasers (comb sources) are not self-referenced to avoid the complexity of amplification and continuum generation. The cw reference lasers are, in our case, stabilized to a high-finesse cavity and calibrated against a self-referenced frequency comb to ˜1 kHz absolute uncertainty. This accuracy carries over to the spectrometer by simultaneous measuring a 1550 nm laser against a self-referenced frequency comb. Such high accuracy matches the resolution and may be required for some very high precision ranging or spectroscopy applications. However, for many coherent LIDAR systems, kHz linewidths at much lower accuracy are adequate and for Doppler-limited spectroscopy ˜MHz accuracy and linewidths are sufficient. Therefore, the high-finesse cavity could be replaced by a simpler cavity or molecular reference line.

FIGS. 10A-10E illustrates a dual-comb spectrometer measurements of a tunable external cavity diode laser (ECDL) swept nonlinearly over a 28 nm span. Unlike approaches with a reference etalon, the output is not continuously tracked; each ˜320 μs measurement is independent of the previous measurements. To obtain the instantaneous frequency in FIG. 10D, we fit the instantaneous phase, θ_(L)(t₂)−δφ₂(t₂), with a 10^(th)-order polynomial to smooth over the ˜0.5 rad per point comb phase noise, and calculate its derivative. As a consequence of limiting our pulse energies to ˜pJ to avoid detector saturation, the comb tooth power is ˜2 nW and the high-resolution spectra's SNR has a shot-noise limit of ˜66 dB. Our SNR is ˜6 dB lower due to added detector noise.

FIG. 11A demonstrates frequency metrology of a mechanically perturbed laser at a time resolution of T_(W)=30 μs.

The coherent dual-comb spectrometer disclosed herein provides for time-resolved measurements of dynamic, arbitrary, cw optical waveforms at ˜320 μs and ˜30 μs time resolution covering a dynamic range in frequency of >10⁹ (the ratio of the 42 nm bandwidth to 3 kHz resolution). The spectrometer can be optimized for time-resolution, optical bandwidth, or precision to match a given application. This technique allows for metrology of dynamic cw lasers and, more generally, circumvents the sensitivity limits of direct comb spectroscopy by combining the precision and accuracy of combs with the high sensitivity and simple detection topology of conventional cw laser sensing.

It should be understood that the operation of the spectrometer can be described in the frequency domain or time domain (see FIGS. 9E-9F and 9A-9C), where it is similar to a linear optical sampling system or Fourier spectroscopy. The operation thereof will first be explained in more detail from a time-domain perspective with a particular emphasis on the wideband spectrum and frequency estimation, and second from a frequency-domain perspective with a particular emphasis on the consequences of heterodyne detection and bookkeeping for the absolute frequency scaling.

In a time-domain picture, the pulse train of comb 1 or 2 is

E _(comb1,2) =e ^(iφ) ^(1,2) Σ_(m) e ^(imθ) ^(1,2) E _(1,2)(t−mT _(1,2)).   (0.1)

where the index 1 or 2 indicates the comb number, θ_(1,2) is the carrier envelope offset (ceo) phase shift, E_(1,2)(t) is the field of a single pulse, φ_(1,2) is an arbitrary phase offset, m is the pulse number, and T_(1,2)=f_(r1,2) ⁻¹ is the repetition period. The difference in pulse periods is ΔT=T₁−T₂, (assumed to be positive). For an input E-field E_(L) (t), the digitized overlap with the m^(th) pulse from comb 1 as a function of sample time t₁=mT₁ is

V ₁(t ₁)=α∫E _(L)(t)e ^(−iφ) ¹ e ^(−it) ¹ ^(f) ^(r1) ^(θ) ¹ E ₁*(t−t ₁)dt,   (0.2)

with an analogous equation for comb 2, assuming a detector bandwidth matched to Nyquist and, for now, IQ detection (so that the complex voltage can be constructed from the sum of the I and Q channels), where a is a constant that converts to ADC counts. We choose the pulse indices such that the pulses from comb 1 and 2 overlap for the same pulse index, m=0. With our phase-locking conditions (see Methods) pulses from the two combs will resynchronize every N comb 1 pulses (which set the digitization clock). We define a single frame as m=0 to N−1, which takes a time Δf_(r) ⁻¹=T_(W), and process on a frame-by-frame basis so that T_(W) sets the time resolution.

Assuming the pulse spectrum varies slowly compared to the laser spectrum, V₁(t₁)=ae^(−iφ) ¹ e^(−it) ¹ ^(θ) ¹ ^(f) ^(r1) {tilde over (E)}₁*(f_(L))E_(L)(t₁) where the tilde represents a Fourier transform. Substituting E_(L)(t)≡E_(L) exp[iθ_(L)(t)] results in

V ₁(t ₁)=aE _(L) exp(i[θ _(L)(t ₁)−φ₁−θ₁ t ₁ f _(r1)])   (0.3)

with a corresponding beat frequency of f₁(t₁)=f_(L)(t₁)−{dot over (φ)}₁(t₁)/(2π)−nf_(r1)−f_(ceo), where the ceo frequency is f_(ceo,1)=(2π)⁻¹θ₁f_(r1).

Returning to Eq. (0.2), and its analog for comb 2, and writing the product V₁(m)V₂*(m) as a function of sample number, m, we find

V ₁(m)V ₂*(m)=a ² ∫∫{tilde over (E)} _(L)(f){tilde over (E)} _(L)*(f′){tilde over (E)} ₁*(f){tilde over (E)} ₂(f′)e ^(i2πmΔTf) ^(opt) e ^(−2πm(fT) ¹ ^(−f′T) ² ⁾ df df′,   (0.4)

(after taking the Fourier transform of each electric field) where the effective offset frequency f_(opt)≡(θ₂−θ₁)(2πΔT)⁻¹ modulo the ambiguity (ΔT)⁻¹ corresponds to the frequency of the tooth of combs 1 and 2 locked to the cw reference laser at 1560 nm (or 1535 nm) for our locking conditions.

To proceed further, we need an expression for {tilde over (E)}_(L)*(f){tilde over (E)}_(L)(f′). For a static cw laser {tilde over (E)}_(L)*(f){tilde over (E)}_(L)(f′)=I_(L)(f)δ(f−f′), where I_(L)(f) is the laser intensity spectrum. The same expression holds true for multiple, uncorrelated cw lasers. An incoherent source will also satisfy this relationship, <{tilde over (E)}_(incoherent)*(f){tilde over (E)}_(incoherent)(f′)>=I_(incoherent)(f)δ(f−f′), after averaging as indicated by the bracket. With this condition, and furthermore assuming the combs have similar spectral filtering and negligible differential chirp (by roughly matched fiber lengths) so that {tilde over (E)}₁(f){tilde over (E)}₂*(f)≡I_(comb)(f), the discrete Fourier Transform (DFT) of V₁(m)V₂*(m) over a single frame is

{tilde over (V ₁ V ₂*)}(n)=a ² I _(comb)(f _(opt) +nf _(r2) +qΔT ⁻¹)I _(L)(f _(opt) +nf _(r2) +qΔT ⁻¹)   (0.5)

where n is a frequency index and runs from 0 to N−1 and q is an integer. (With our phase locking conditions, Nf_(r2)=ΔT⁻¹). The above would be true if we had full IQ detection; instead we have heterodyne detection and Eq. (0.5) is modified to equal {tilde over (V)}{tilde over (V₁V₂*)}(n)=a²I_(comb)(f_(opt)±nf_(r2)+qΔT⁻¹)I_(L)(f_(opt)±nf_(r2)+qΔT⁻¹) with n running from 0 N/2. The DFT of V₁(m)V_(2S)*(m), where V_(2S)(m) are the samples generated from the shifted comb 2, is {tilde over (V)}{tilde over (V₁V_(2S)*)}(n)=a²I_(comb)(f_(opt)+(4ΔT)⁻¹±nf_(r2)+qΔT⁻¹)I_(L)(f_(opt)+(4ΔT)⁻¹±nf_(r2)+qΔT⁻¹) since a frequency shift of 1.25 f_(r2) of comb 2 effectively increases f_(opt) by (4ΔT)⁻¹.

The expression {tilde over (E)}_(L)*(f){tilde over (E)}_(L)(f′)≈I_(L)(f)δ(f−f′) is not true for a frequency-modulated cw laser and spurious tones will appear in the calculated intensity spectrum similarly to conventional Fourier spectroscopy. The spurious tones are suppressed with averaging of multiple data frames, provided the frame rate is not phase coherent with the modulation. More relevant, the above expression is also not true for a strongly swept cw laser, i.e. one that changes by >f_(r) over the interval T_(W). We discuss this effect and its consequences next.

With the assumption that T_(W) is reasonably short, we approximate any laser frequency change as linear, so that

θ_(L)(t)=θ₀+2πf _(L0) t+πCt ²   (0.6)

in terms of a linear chirp C and starting laser frequency f_(L0). The laser frequency is f_(L)(t)≡{dot over (θ)}_(L)(t)/(2π)=f_(L0)+Ct so that it covers a bandwidth CT_(W) over the interval T_(W). From Eq. (0.3) and assuming zero ceo phase as in the main text and zero carrier phase noise, the product V₁(t₁)V₂*(t₂)∝|E_(L)|² exp[i(θ_(L)(t₁)−θ_(L)(t₂))]=|E_(L)|² exp[i(2πf_(L0)Δt+πC(t₁ ²−t₂ ²))], where Δt=t₁−t₂. We can rewrite this equation as V₁(t₁)V₂*(t₂)∝|E_(L)|² exp(iθ′_(L)(Δt)), where the phase factor is θ′_(L)(Δt)≈2πf_(L0)Δt+2πNCΔt² and N=t₂/Δt=T₂/ΔT=f_(r1)/Δf_(r) is the down-sampling factor (see Methods). Substituting Δt→t and NC→C to account for the down-sampled time, the phase θ′_(L) (Δt) equals Eq. (0.6), except that there is an additional factor of 2 in the last term describing the chirp. Therefore, the calculated wideband spectrum (i.e. spectrum of V₁(t₁)V₂*(t₂)) over a frame of length T_(W) will have a width ˜2 CT_(W), twice as wide as it should be. More significantly this increased width all occurs at positive frequencies (for an increasing chirp) which leads to a bias in the calculated mean frequency equal to CT_(W)/2. We have verified this effect by plotting the calculated wideband spectrum for a continuously swept laser; we find an overlap of consecutive wideband spectra and a spectral width roughly twice the expected value given the chirp rate.

For chirp rates C<f_(r)/T_(W), the shift is less than half the wideband spectrum resolution of f_(r), and can be ignored. However, if a chirp C>f_(r)/T_(W), we correct the estimated frequency by CT_(W)/2, where we calculate the chirp C directly from the individual frequencies f_(1,2).

For most situations, this corrected mean frequency is more than adequate. However, it would, in principle, be possible to track the instantaneous frequency over all time even for high chirp without correction by simply successively switch between f₁, f₂, and f_(2S) as the laser frequency crosses successive comb teeth of comb 1, 2 and the shifted 2.

The operation can also be described in the Frequency domain. We are interested here in the frequency assignment of the high-resolution spectrum. In the frequency comb description, comb 1 is E_(comb,1)(t)=Σ_(n)|{tilde over (E)}₁(f_(1n))|e^(−i(2πf) ^(1n) ^(t+φ) ^(1n) ⁾ where the nth comb tooth has phase φ_(1n) and frequency f_(1n)=f_(opt)+nf_(r1) as measured relative to the frequency of a specific, phase-locked tooth at f_(opt) (rather than the ceo frequency). Similarly, the teeth of comb 2, labeled by index n′, have frequency f_(2n′)=f_(opt)+n′f_(r2) and phase φ_(2n′) (with respect to the same common phase-locked tooth at f_(opt) for our phase-locking conditions). The heterodyne signal between the laser and the nearest tooth of combs 1 and 2, digitized synchronous with comb 1 (i.e. at a rate f_(r1)), is

V ₁(t)=√{square root over (2)}a|{tilde over (E)} ₁(∫_(1n))∥E _(L)(t)|cos(±2π∫₁ t)

V ₂(t)=√{square root over (2)}a|{tilde over (E)} ₂(f_(2n′))∥E _(L)(t)|cos(±2π(f ₂ +zf _(r1))t)   (0.7)

where the arguments of the cosine function are 2πf₁t=θ_(L)(t)−2πf_(1n)t−φ_(1n) and 2πf₂t=θ_(L)(t)−2πf_(2n′)t−φ_(2n′) and we have let t₁→t to minimize subscripts. As a consequence of the experimental reality that the signal from the second comb is sampled at t≡t₁ rather than t₂, aliasing can take place for f₂, which is accounted for by zε{−1,0,1}, as needed to force |f₂+zf_(r1)|≦f_(r1)/2.

Without IQ detection, the corresponding “sine” component is unavailable and the product V₁(t)V₂(t) has both a difference and sum frequency term. Substitution of f₁, f₂ and n=n′−z (i.e. assuming the common lock point f_(opt) is near f_(L)), the difference frequency term is [V₁(t)V₂(t)]_(diff)=a²I_(comb)(f_(2n))I_(L)(f_(L))cos[2πnΔf_(r)t−Δφ₁₂] where Δφ₁₂=φ_(1n)−φ_(2n) and letting n′→n. In other words, the amplitude at the rf frequency nΔf_(r) is proportional to the input intensity at f_(1n)≈f_(2n), or more generally, the intensity within the spectral band f_(2n)±f_(r2)/2 closest to the n^(th) tooth of comb 2. When scaled up in frequency by the down-sampling factor N, the rf spectrum is the wideband intensity spectrum. (This scaling is automatic if we rescale the spacing between time points to ΔT before the Fourier transform.) A similar signal is generated using comb 1 and the shifted comb 2, after shifting the comb 2 teeth by f_(2n)→f_(2n)+f_(r2)/4. Note that residual phase noise between the combs associated with Δφ₁₂ must be minimized to avoid broadening of the difference frequency signal beyond Δf_(r), with a corresponding loss of resolution.

The sum frequency gives a very different, and undesired, result of [V₁(t)V₂(t)]_(sum)=a²I_(comb)(f_(2n))I_(L)(f_(L))cos[2θ_(L)(t)−4πf_(opt)t−4πf_(r2)t−φ_(1n)−φ_(2n)], which depends strongly on the relative phase between the input laser and comb (and therefore reference laser). For differential phase noise >>1 rad (i.e. for our narrow comb lines a laser linewidth greater than the time-bandwidth product) this sum frequency signal will be broad and easily distinguishable from the difference frequency signal, which we can then use to assign an absolute frequency to the high-resolution spectrum. If the instantaneous laser linewidth is below f_(r)/2, a Hilbert transform can also be applied to eliminate the undesired sum frequency contribution. For an instantaneous laser linewidth greater than f_(r)/4, the sum frequency signal is so broad as to be distributed evenly from 0 to f_(r)/2 and effectively becomes the limiting noise source. Note that for very short update intervals, it might be that θ_(L)(t)−2πf_(L)t does not vary more than ˜1 rad and the sum-frequency signal would be narrow and therefore indistinguishable from the difference-frequency signal (although that case was never met here even for a DFB fiber laser at T_(W)=300 μs). In that case, a common-mode phase modulation can be applied to the combs to purposefully broaden the sum frequency, while maintaining negligible residual phase noise between the combs (i.e. Δφ₁₂=0).

It should be understood that like reference numerals identify corresponding or similar elements throughout the several drawings. It should also be understood that although a particular component arrangement is disclosed in the illustrated embodiment, other arrangements will benefit herefrom.

Although particular step sequences are shown, described, and claimed, it should be understood that steps may be performed in any order, separated or combined unless otherwise indicated and will still benefit from the present disclosure.

The foregoing description is exemplary rather than defined by the limitations within. Various non-limiting embodiments are disclosed herein, however, one of ordinary skill in the art would recognize that various modifications and variations in light of the above teachings will fall within the scope of the appended claims. It is therefore to be understood that within the scope of the appended claims, the disclosure may be practiced other than as specifically described. For that reason the appended claims should be studied to determine true scope and content. 

1. A method of comb-based spectroscopy with synchronous sampling for real-time averaging comprising: generating a source comb that transmits a source pulse train through a sample to provide a signal pulse train; generating a local oscillator (LO) comb that transmits a LO pulse train at a different repetition rate than the source pulse train so that the difference in repetition rates between the combs is Δf_(r); establishing a base coherence between the source comb and the LO comb wherein each Nth pulse from the LO pulse train will have identical overlap with each Nth+1 pulse from the source pulse train in terms of both the pulse envelope and the carrier phase; detecting the signal pulse train through linear optical sampling against the LO pulse train in which consecutive samples of an overlap between a signal pulse and a LO pulse from the LO pulse train yields a measurement of the signal pulse; digitizing the overlap between the signal pulse train and the LO pulse train synchronously with the LO pulse train; and real time summing of every Nth digitized sample to generate an averaged signal pulse.
 2. A method as in claim 1, wherein the coherence between the LO comb and source comb is established by phase locking a first tooth of the source comb and LO comb to a first CW reference laser and a second tooth of the source comb and the LO comb to a second CW reference laser with the number of intervening teeth for the source comb being N1 and the number of intervening teeth for the LO comb being N2 such that N2/(N1−N2) is an integer N.
 3. A method as recited in claim 1, further comprising: optically filtering the source and LO comb outputs so that Δf_(r) can be high enough to avoid 1/f noise without violating the Nyquist sampling constraints; successively shifting the center filter position across the full LO and source comb spectra to measure the signal pulse in different spectral regions; coherently reconstructing the signal pulse by concatenating the measurement at different filter positions in the frequency domain; shifting one of the source comb and the local oscillator (LO) comb at each filter position by a multiple of Δf_(r) such that the carrier frequency of the averaged signal pulse is moved away from base band or Nyquist.
 4. A method as recited in claim 1, wherein: a portion of the source pulse train that does not pass through the sample provides a reference pulse train; delaying the reference pulse train relative the signal pulse train by a percentage of a period of the signal pulse train such that the reference pulse train and the signal pulse train do not overlap to provide a resultant interleaved signal pulse train and reference pulse train; and detecting the resultant interleaved signal pulse train and reference pulse train through linear optical sampling against the LO pulse train as in claim 1 to generate a measurement of the interleaved signal and reference pulses.
 5. A method as recited in claim 4, wherein: separating the resultant measurement of the interleaved signal pulse and reference pulse to normalize the signal pulse by the reference pulse yielding the time-domain response of the sample.
 6. A method as recited in claim 5, wherein the Fourier transform of the normalized signal pulse yields the complex frequency domain response of the sample.
 7. A method as recited in claim 1, wherein: a portion of the source pulse train that does not pass through the sample provides a reference pulse train; and alternately detecting the signal pulse train and the reference pulse train through selectively blocking each of the signal pulse train and the reference pulse train to normalize the signal pulse by the reference pulse yielding the time-domain response of the sample.
 8. A method of comb-based spectroscopy for measuring a CW source at time-bandwidth limited resolution by using frequency combs with a high degree of mutual coherence (<1 radian phase noise) comprising: generating a first comb that transmits a first pulse train that is optically combined with a CW source; detecting the overlap of each pulse of the first comb with the CW source in a first photodetector; digitizing a photodetector response from the first photodetector for each pulse of the first pulse train; generating a second comb that is optically coherent with the first comb, said second comb transmits a second pulse train that has a pulse period that differs by ΔT and a pulse repetition frequency that differs by Δf_(r) from the first pulse train, and is separately combined with the CW source; detecting the overlap of each pulse of the second comb with the CW source in a second photodetector; digitizing a photodetector response from the second photodetector for each pulse of the second pulse train; multiplying the digitized samples from the first comb times the digitized samples of the second comb over a time (1/Δf_(r)) to generate a data point record of length (f_(r)/Δf_(r)) where f_(r) is the repetition rate; rescaling the point separation for the data record therebetween to be ΔT; and Fourier transforming the product of the digitized samples to yield a wideband spectrum of the CW source at a resolution given by the comb repetition rate and an ambiguity given by 1/(2ΔT) to provide an absolute frequency of the CW source with respect to a CW reference laser to which the respective first comb and the second are locked.
 9. A method as recited in claim 8, further comprising: generating a third comb by splitting comb 2 into two parts and frequency-shifting one part to create a third comb with the same repetitition rate but an offset frequency shift equal to a fraction, X, of the repetition rate; detecting the overlap of each pulse of this third frequency-shifted comb with the CW source in a third photodetector; processing the digitized samples from the third comb and the digitized samples from the first comb as described in claim 8 to yield a wideband spectrum of the CW source at a resolution given by the comb repetition rate and but an offset frequency that differs by X/ΔT from the wideband spectrum generated from the product of the first and second combs in claim
 8. Combining the two wideband spectrum to provide an absolute frequency of the CW source with respect to a CW reference laser to which the respective first comb and the second are locked to an overall ambiguity 1/(ΔT) and without ambiguities related to the Nyquist frequency 1/(2ΔT)
 10. A method as recited in claim 8, further comprising: individually Fourier transforming the digitized samples from either of the first comb or the second comb or the third comb to generate a spectrum of the CW source with time-bandwidth limited resolution.
 11. A method as recited in claim 8, wherein the phase coherence between the source and LO combs is established as recited in claim
 2. 